Prom incommensurate correlations to mesoscopic spin resonance in YbRh 2 Si2 
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Spin fluctuations are reported near the magnetic field driven quantum critical point in YbRli2Si2. 
On cooling, ferromagnetic fluctuations evolve into incommensurate correlations located at qo=± 
(5,8) with 5=0.14 ± 0.04 r.l.u. At low temperatures, an in plane magnetic field induces a sharp 
intra doublet resonant excitation at an energy Eo = gps^oH with g=3.8±0.2. The intensity is 
localized at the zone center indicating precession of spin density extending £ = 6 ± 2 A beyond the 
4/ site. 

PACS numbers: PACS numbers: 75.40.Gb, 74.70.Tx, 75.50.Cc 
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The development of inter site coherence amongst 
Kondo-screened ions is a phenomenon that is central 
to understanding the heavy fermion state and quantum 
phase transitions. [l|-|4j Here we use neutron scattering to 
monitor such correlations as a function of magnetic field 
and temperature in a clean system where disorder due 
to chemical substitution is not expected. These exper- 
iments show the formation of a coherent Kondo lattice 
and a low temperature field driven transition to a dressed 
single ion phase. 

YbRh 2 Si2 is a heavy fermion metal with weak antifer- 
romagnetic order at very low temperatures ~ 50 mK. [5] 
The resistivity varies linearly with temperature (p ~ T) 
demonstrating strong non- Fermi liquid character. Q In 
a small magnetic field of ~ 0.7 T along the c axis (and 
^0.06 T when the field is applied within the a — b plane), 
the magnetic order is suppressed and the non-Fermi liq- 
uid phase is transformed continuously into a Landau 
Fermi liquid where p ~ T 2 . 0, [j| The large change in Hall 
number and band structure calculations indicate a quan- 
tum critical transition in the electron density which can 
be interpreted in terms of Fermi surface reconstruction or 
local quantum criticality. UlSl At higher fields de Haas 
van-Alphen measurements indicate a Lifshitz transition 
of the Fermi surface. 
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Despite the heavy fermion nature, electron-spin reso- 
nance (ESR) shows a signal indicative of localized Yb 3+ 
4/ moment behavior under the application of a field. (IBj 
While some work has been performed on the magnetic 
dynamics, little is known about the magnetic response 
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across the field tuned quantum critical point. 

We present a neutron inelastic scattering study of sin- 
gle crystalline YbRli2Si2. We will show that at zero field, 
the magnetic fluctuations are incommensurate. Under 
the application of a magnetic field, driving the system 



through the quantum critical point, commensurate un- 
derdamped fluctuations develop at the Zeeman energy of 
the Yb 3+ 4/ crystal field doublet. The momentum lo- 
calized resonance, however, indicates a spin density that 
extends beyond the 4/ ion. 

Experiments on single crystalline YbRh.2Si2 prepared 
using Zinc flux were performed on the SPINS and OSIRIS 
spectrometers and powder experiments (discussed in the 
supplementary information) used MARL Approximately 
300 (with total mass ~ 3 g) single crystal samples (growth 
and characterization described in Ref. 1181 ) were coaligned 
in the (HHL) scattering plane using a similar method as 
described previously. [19j Experiments on SPINS used a 
focussing analyzer with 11° acceptance, =3.7 meV and 
vertical magnetic fields up to 11 T. To probe the low en- 
ergy dynamics (E <0.5 meV) we made use of the OSIRIS 
indirect spectrometer with E/=1.84 meV. The data on 
OSIRIS has been corrected for a background obtained 
by imposing detailed balance. [20j j In both experiments 
the temperature was monitored by a resistive thermome- 
ter attached near the sample. The lowest temperature 
achieved was 100 mK, therefore this experiment does not 
access the ordered state. Absolute normalization was 
performed against the (0,0,2) nuclear Bragg peak and 
through a comparison of the reported magnetic field de- 
pendent magnetization. [2l| 

We first describe the magnetic fluctuations at low tem- 
peratures near the ordering transition. Figure Q] shows 
constant energy scans near the ferromagnetic Q=(0,0,2) 
position at 100 mK. Panel (a) summarizes the OSIRIS 
and SPINS data showing ridges extending along Two. This 
demonstrates that the momentum dependence is con- 
trolled by a higher energy scale such as the Kondo tem- 
perature or the Fermi energy. As shown in the supple- 
mentary information, these incommensurate low-energy 
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FIG. 1. The Tiuj — Q dependence of the magnetic fluctua- 
tions in YbRh2Si2 at T=0.1 K. (a) the position in Q of the 
ridges of scattering where filled (open) circles are from SPINS 
(OSIRIS). (6) and (c) plot contour maps of the spin fluctua- 
tions taken on OSIRIS at 0.3 K and 15 K integrating over the 
range (q,q,2±0.25). (d) and (e) show background corrected 
cuts in momentum at 0.5 meV. The horizontal bars are the 
experimental resolution. 



fluctuations are well below the first crystal field level lo- 
cated at 17.9 ± 0.6 meV and therefore are associated with 
allowed transitions within the ground state doublet. 

Panels (d) and (e) show representative scans taken 
on SPINS at 0.5 meV with a constant background sub- 
tracted. The scan along the (H,H,0) direction (d) illus- 
trates an incommensurate modulation in the basal plane 
while the scan along L (panel (e)) shows ferromagnetic 
inter plane correlations. An (H,H,2±0.25)-?iw slice (panel 
(b)) demonstrates that the incommensurate ridges extend 
to the lowest energy transfers accessed (~ 0.1 meV). The 
fluctuations are peaked at = (5,5) with 5 =0.14 ± 
0.4 with no measurable dispersion or offset along L. We 
therefore expect (5, 5) to be the characteristic wavevec- 
tor of the low temperature spin density wave order corre- 
sponding to the transition reported using [j,SR and sus- 
ceptibility. 0, 0] While there is no conclusive evidence 
that Fermi surface nesting drives this transition, the ob- 
served critical wavevector is not far from the spacing be- 
tween large areas of the computed Fermi surface. [l2[ Low 
temperature diffraction and calculations of x(Q) wn l be 
required for progress. 
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FIG. 2. (a) illustrates the linewidth of the magnetic excita- 
tions taken from constant Q scans represented in panels (6) 
and (c). The dashed curve in (b) centered at E=0, is the mea- 
sured resolution on SPINS. Panel (d) summarizes all of the 
temperature dependence from SPINS and OSIRIS and plots 
fcsT^x" as a function of E/ksT over the temperature range 
T=[0.1, 30] K (as described by the colorbar) and integrating 
over Q=(±0.1, ±0.1, 2±0.25) on OSIRIS. The exponent a was 
fitted to be a=1.05 ± 0.03 as described in the supplementary 
information. 



At elevated temperatures of (15 K, panel (c)), the in 
plane response changes considerably with the incommen- 
surate fluctuations being replaced by commensurate fer- 
romagnetic excitations forming a ridge at the zone cen- 
ter. These results confirm a competition between ferro- 
magnetic q=0 excitations and incommensurate spin fluc- 
tuations as inferred from NMR based on a c omp arison 
between the Knight shift and relaxation rate. [22j 

A more detailed survey of the T dependent spin fluc- 
tuations can be found in Fig. [2] Constant Q=(0,0,2) 
scans at 30 and 6.5 K are shown in panels (b) and 
(c). The solid curves are a fit to the relaxational form 
X"(Q,W) = XQ^r(T)/(r(T) 2 + c 2 ), where T(T) is the 
relaxation rate and x'q is the susceptibility. The tem- 
perature dependence of T(T) (Fig. [5] (a)) varies linearly 
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FIG. 3. (a) illustrates a background subtracted constant Q 
scan showing the presence of a sharp resonance peak on the 
application of a 5T field on SPINS, (b) and (c) illustrate con- 
stant energy scans at the resonance energy and below at 0.5 
meV. (d) illustrates a contour taken on OSIRIS illustrating 
that the resonance peak at 1.25 T is sharp in H and energy. 



and reaches zero close to T=0. Such behavior is expected 
near quantum criticality. 

In a similar manner to the case of scaling near a clas- 
sical phase transition driven by thermal fluctuations, we 
investigate the scaling properties of the dynamic suscep- 
tibility near the low- field quantum critical point. Several 
scaling theories have been proposed for the susceptibil- 
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ity. [13j, |23j The uniform dynamic susceptibility \" a t 



all temperatures and energies from SPINS and OSIRIS 
are compiled in (d) which shows (fcsT) a x" as a function 
of E/k B T. The data is best fit with a =1.05 ± 0.03 
consistent with 1 as expected in the case of an itinerant 
ferromagnetic at T>T 'curie and T oc 1/Xq, when x'q na s a 
Curie form. Q The exponent is the same as observed for 
the scaling of heat capacity. [7| Such a scenario is consis- 
tent with the fluctuations being close to ferromagnetic as 
reported in Fig. [T]and are similar to the case of UCoGc 
which is close to a critical point between ferromagnetism 
and superconductivity. [24[ 

We now investigate the possibility of anomalous scaling 
in our experiment which would manifest as a deviation 



FIG. 4. (a) illustrates the magnetic field dependence of the 
resonance energy and linewidth. (b) shows the integrated in- 
tensity in absolute units. Representative constant Q scans are 
shown in panels (c)-(e) at various magnetic fields. 



from the conventional scaling described above. [23j A fit 
to k B T^ X " a s a function of E/(k B T) resulted in /3=0.92 
±0.05 also consistent with 1 and different from other sce- 
narios used to describe the fluctuations in CeCu2Si2 and 
CeCue-^Auz where such an analysis gave /3=1.5. [25| 
The critical dynamics in CeCug _xAu a was also found to 
be described with cv=0.75 (Ref. I26T I. also different from 
YbRli2Si2- While anomalous scaling may exist close to 
the quantum critical point, our data at Q ~ does not re- 
quire such forms. Indeed, the scaling exponent we derive 
indicates that we are in the quantum disordered regime 
and higher resolution measurements at fields, tempera- 
tures, and wavevectors closer to the critical point are 
required to make contact with the anomalous thermody- 
namic data. 

Next we study the low temperature incommensurate 
fluctuations as a function of field. The result is summa- 
rized in Fig. |31 for an [110] field where the critical field 
is ~ 0.06 T. Fig. [3] (a) illustrates that at 5 T a spin 
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resonance develops with the solid line a fit to an under- 
damped simple harmonic oscillator plus a relaxational 
form to describe the low-energy dynamics, (b) and (c) 
plot constant energy scans at the resonance energy (1.05 
meV) and below at 0.5 meV. The resonance is found to 
be localized in Q space with the fluctuations below the 
resonance retaining the incommensurate structure mea- 
sured at zero field. The solid curve in (b) is a resolution 
convolved fit to a lorentzian squared resulting in a dy- 
namic correlation length of £=6 ± 2 A. This indicates 
that the resonance involves a finite region of correlated 
spins, (d) illustrates a contour plot taken on OSIRIS at 
1.25 T demonstrating that the resonance peak (even at 
small magnetic fields) is localized in energy and momen- 
tum. Searches for dispersing spin excitations located in 
q away from the resonance peak, as measured in ferro- 
magnetic MnSi (Ref. [27]) and calcium ruthenate (Ref. 
28T ) under an applied magnetic field in the paramagnetic 
state, failed to observe any dispersing modes. There- 
fore, the resonance peak measured here is not indicative 
of field induced spin-waves, like in other itinerant ferro- 
magnets, but rather represents a coherent precession of 
spin density extending beyond a 4/ site. 

The magnetic field dependence of the resonance peak is 
displayed in Fig. |4j (a) shows that the resonance energy 
varies as gfj,B^oH with g=3.8 ±0.2 and remains under- 
damped at all measured fields with little change in the 
linewidth. Background corrected constant Q scans track- 
ing the resonance with field are shown in panels (c)-(e). 
The integrated weight is also field independent as shown 
in (b) and matches the calculated spectral weight for intra 
doublet transitions (or transverse fluctuations) presented 
in the supplementary information. Therefore all of the 
Yb 3+ ions are contributing. The g factor is comparable 
to 3.6 obtained from ESR (Ref. [3) and is consistent with 
the crystal field scheme described in the supplementary 
information. Therefore, the ESR signal and the neutron 
scattering resonance have a common origin and are both 
associated with the field split ground state doublet of the 
Yb 3+ ions. 

The presence of both a sharp magnetic field induced 
resonance and an ESR signal are unique to YbRli2Si2 
amongst heavy fermion materials. In heavy electron sys- 
tems, an ESR signal is typically only observed in the 
presence of local Kondo impurities. At the phemenolog- 
ical level, the neutron experimental result is similar to 
the Haldane spin chain Y2BaNii_ 2; Mg 2 ;05 with Mg im- 
purities. [29| In that system, a magnetic field was found 
to induce a resonance peak which is sharp in energy and 
momentum representing a staggered magnetization near 
the edge of a Haldane chain segment. The finite linewidth 
in momentum was associated with a dynamic correlation 
length which measures the extent of the impurity edge 
state. Similar arguments have been applied to low energy 
resonant field induced effects in underdoped cuprates, 
possibly the result of exciting free spins near charge rich 



regions. [30( 

A similar physical picture maybe applied to YbRh.2Si2. 
A magnetic field induces localized droplets of Yb 3+ spins 
which can be resonantly excited through intra doublet 
transitions. The localized region of spins are analo- 
gous to a Kondo impurity which would give rise to an 
ESR signal. The length scale for such a Kondo impu- 
rity spin, should be related to £ ~ Kvf /ksl 'Kondo = 
h 2 k f /mk B T Kondo . [3l| Taking TV orldo =24 K, k f -0.5 
a*, and 7 - 1.5 J/mol K 2 (Ref. at 0.1 K at zero field) 
yields a lengthscale of — 15 A. The result is on the same 
order as the dynamic correlation length of 6 A extracted 
above indicating that £ maybe set by the Kondo temper- 
ature and Fermi velocity. Because the effect is not due 
to purely localized Yb 3+ ions and is not associated with 
correlated dynamics over long length scales, we refer to 
this resonance as mesoscopic. The lack of Fermi surface 
nesting at high fields, the large Fermi velocity, and the 
heavy nature of the bands appears to result in the for- 
mation of Kondo screened ions. Such effects appear to 
be supported by ferromagnetic fluctuations and similar 
results may exist in other ferromagnetic heavy fermion 
metals. [jillallH! 

Our experiments indicate that the low temperature 
zero field dynamics of YbRli2Si2 is effected by Fermi sur- 
face nesting while at high fields and temperatures the 
response mimics Kondo screened spins above a ferromag- 
netic transition. 
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Supplementary information for "From incommensurate correlations to mesoscopic 

spin resonance in YbRh 2 Si2" 

(Dated: July 19, 2012) 

We present supplementary information regarding the sample preparation and experimental details, 
crystal field analysis used to provide a benchmark for the neutron scattering cross sections, and 
further information regarding the scaling analysis presented in the main text. 
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Supplemental material describing experimental 
setup and samples 

One of the largest hurdles to obtaining neutron inelas- 
tic scattering data on many heavy fermion systems is the 
possibility to synthesize large single crystals on the or- 
der of several grams. We have overcome this obstacle 
in YbRli2Si2 by coaligning ~ 300 samples on aluminium 
plates using the facets and the plate like geometry of the 
single crystals. Each sample was fixed to the plates with 
Fomblin hydrogen-free oil and further secured by wrap- 
ping the plate in aluminum foil. The resulting rocking 
curve and sample mount are illustrated in Fig. EU The 
rocking curve was measured on SPINS using a configu- 
ration of guide - PG(004) - 80 - S - 80 - detector at 
room temperature with Ej=14.8 meV. 

The individual samples were prepared by the same 
method as described in Ref. QJ The samples used here 
were grown from Zinc flux which facilitated the growth of 
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FIG. SI. The rocking curve measured at Q=(0,0,4) taken 
on SPINS with Ei = 14.8meV. The inset illustrates the Alu- 
minium foil covered plate geometry of the sample mount. The 
c-axis was aligned perpendicular to the plate face. The crys- 
tals were grown using a Zinc flux as described in Ref. |l|. 



larger samples than produced using indium. A detailed 
comparison of the magnetic properties, between samples 
prepared using both both methods, characterized with 
ESR is in Ref. 0. It was shown that, while the resid- 
ual resistivities of the zinc prepared samples are larger 
than in materials used in early studies, indium flux typi- 
cally produces higher resistivity ratios and sharper ESR 
linewidths. Heat capacity studies (Ref. 1) did show that 
the samples prepared using zinc displayed the same diver- 
gence as reported for indium prepared samples. The sus- 
ceptibility at temperatures above ~ 1 K was also found 
to be in good agreement between samples prepared using 
zinc and indium. 



Supplemental material describing crystal fields 

In this supplemental section, we present a self consis- 
tent description of the crystal field energies, intensities, 
Lande g factors along with the small ordered moment and 
the integrated transverse spectral weight reported above. 
A crystal field analysis has been discussed in Refs. and 
4, however, for completeness, a similar analysis is shown 
here in the context of our neutron measurements with 
the goal of understanding the integrated spectral weight. 

The crystal field excitations were measured on the 
MARI chopper spectrometer with incident energies of 50 
meV and 100 meV. The powder sample was loaded into a 
thin annulus to minimize effects due to absorption. Fur- 
ther scans using Ei=200 meV found no scattering at en- 
ergies above 50 meV energy transfer. To model the back- 
ground from phonons, we have performed similar scans 
at large wave vector transfers (integrating over [10, 12] 
A" 1 and [6.5, 7.5] A^ 1 for E;=100 and 50 meV respec- 
tively) at 150 K where the scattering is dominated by 
lattice vibrations and the Yb 3+ form factor ensures min- 
imal magnetic scattering. These scans were then scaled 
to the low temperature and low-Q data such that low 
and high energy transfers were subtracted to zero. The 
resulting data in Fig. [S2] illustrating three distinct crys- 
tal field excitations will be the basis for the subsequent 
analysis. 

For Yb 3+ with J = 7/2 in a local environment, 
the crystal field scheme should consist of 4 doublets. The 
corresponding Hamiltonian takes the following form, 
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FIG. S2. The crystal field scheme measured on MARI at 5 K. 
(a) is a schematic illustration of the ground state consisting of 
4 doublets, (b) and (c) show constant Q cuts from MARI data 
at Ei=100 meV (50 meV), integrating over Q=[1.5, 2.5] A -1 
([2, 3] A -1 ). The data have been corrected for a background 
resulting from phonons taken at high momentum transfers 
and temperature. 



where are the Stevens operators. [3] 

To calculate the experimental values of the Stevens 
parameters (B™), we have fit the data in Fig. [S2] to 
three Gaussians to extract three integrated intensities 
and three energies. The crystal field excitations are 
broader than the resolution of the spectrometer - partic- 
ularly at higher energy transfers where the resolution on 
MARI improves over lower energy transfers. This might 
be indicative of a finite lifetime for the excitations, or pos- 
sibly the considerable dispersion noted in recent ARPES 
studies. @ For the lower two crystal field levels, we have 
averaged the results from the Ei=100 and 50 meV data 
sets. An additional data point from neutron scattering is 
provided by Fig. 4 from which an in-plane g factor was 
extracted. We have also included a g factor parallel to 
c described in Ref. This resulted in 8 "data points" 
from which the 5 Steven's parameters were fit. The result 
is tabulated below. 

An analytic analysis found that the energy eigenvalues 
of the crystal field Hamiltonian only depends on cross 
terms proportional to the product of B\ and B|. This 
introduces a sign ambiguity of both parameters in our 
fitting analysis and is emphasized in the table by (±). A 




FIG. S3, a) — c) illustrate a series of scaling plots for vari- 
ous fixed values of a over all temperatures and energy range 
studied, d) shows a plot of \ 2 as a function of a demonstrat- 
ing that the minimum (denoted by the dashed line) exists at 
a=1.05. The concavity is a measure of the errorbar. 



Bi 
B° 6 
Bi 



0.58 ± 0.10 meV 
0.0122 ± 0.008 meV 
(±) -0.047 ± 0.011 meV 
-3.0 e-4± 0.8 e-4 meV 
(±) 0.0102 ± 0.0016 meV 



point charge (Ref. |6|) analysis summing over a cluster of 
~ 1000 unit cells shows that B\ has the opposite sign to 
that of -Bi; and confirms that B\ and B\ have opposite 
signs. While the applicability of the point charge model 
to YbRli2Si2 is debatable owing to strong bonding (Ref. 
it indicates B\ is negative and B\ is positive. 
Based upon the parameters listed above, the calculated 
positions of the three doublet excitations are 17.9, 25.2, 
40.8 meV and the in-plane g factor is 3.4. These should 
be compared with the measured values of 17.9 ± 0.6, 25.2 
± 0.6, 40.8 ± 0.9 meV and the in plane g factor of 3.8 ± 
0.2. The g factor along c is predicted to be 0.20 consis- 
tent with 0.17 ± 0.07 reported by electron spin resonance 
measurements. Q The resulting ground state doublet has 



3 



the following form (for the two solutions listed above), 

|0+> = (-0.77 ± 0.05)||> + (-0.63 ± 0.05)| - |), 
|0_) - (0.63 ± 0.05)||) + (0.77 ± 0.05)| - |>, 

This results in a low temperature ordered moment of 
<7|{0| J z \0}\fJ.B — 0.010 ± 0.003/i_B and a transverse cross 
section (corresponding to intra doublet transitions) of 
T, l=x , y IH^I+MVs = 5-8 ± 1.0/1%. The small ordered 
moment is consistent with [iSR measurements suggest- 
ing values of ~ 10~ 2 -10~ 4 [Ib- [13] The transverse cross 
section agrees within error with the value reported above 
to reside in the field induced resonance peak. This cor- 
roborates our claim that the field induced resonance re- 
sults from transverse spin fluctuations within the lowest 
energy doublet. 



CeCu2Si2 and CeCue-zAuz) and therefore it is impor- 
tant to quantify how well our data fits to the conventional 
scaling formalism described above. Such deviations are 
quantified in the exponent a and Fig. IS3I illustrates how 
the errorbar on a was derived. Panels a — c) show how 
the data scale as a function E/ksT for several different 
exponents a. Panel c) demonstrates that while the high 
temperature data are very well described with a=0.75, 
the low temperature data, closest to the critical point, 
deviate significantly from this scaling form. A plot of 
the normalized \ 2 ; a measure of the quality of the scal- 
ing analysis, demonstrates that a=1.05 provides the best 
fit to the data. From the concavity, we derive an error- 
bar ±0.03. Therefore, based on this analysis we do not 
observe anomalous scaling in the magnetic dynamics in 
YbRli2Si2 and find the temperature dependence to be de- 
scribed by relaxational dynamics expected in an itinerant 
fcrromagnet. 



Supplemental material describing scaling analysis 

In this section we provide supplementary details of the 
scaling analysis described above. The scaling analysis 
is presented over a broad range in temperature (0.1 K- 
30 K) with the data from SPINS derived from constant 
Q scans at Q=(0,0,2) and the data from OSIRIS inte- 
grated over Q=(±0.1,±0.1,2±0.25). The experiments 
were performed at temperatures above the magnetic or- 
dering transition expected at ~ 0.05 K. No magnetic 
Bragg peak or diverging correlation length was observed 
over this temperature range and therefore it is not clear 
whether the experiment is in the critical regime to ob- 
serve anomalous critical dynamics. We do find the 
temperature dependence of the magnetic dynamics to 
be well represented by a relaxational form x"(Q, w) = 
x' Q wr(T)/(r(r) 2 +w 2 ), where T(T) is the relaxation rate 
and Xq is the susceptibility. For an itinerant ferromag- 
netic at T>T ' curie, T oc 1/x'q an d Xq nas a Curie form 
(oc 1/T). The fact that all of the data collapse onto this 
curve implies that this scaling form is correct. 

Deviations from this conventional scaling form have 
been reported in other heavy fermion systems (namely 
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